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C\l ■ Abstract 

In rotationally constrained percolation models, a site of a percolation cluster could be occupied 

l/^ \ more than once from different directions due to the nature of the rotational constraint. A state 

variable Si is assigned to each lattice site whose value corresponds to the number times it has been 

^^ . visited during the growth of a cluster. It is proposed here that the percolation transition and the 

multifractal aspects of infinite percolation clusters under rotational constraint can be studied defin- 
ed 
C^ ' ing suitable measures in terms of the state variable Sj. This method does not require to introduce 

n^ • any external agency like an electric current or a random walker in order to explore multifractality as 

a, 

O . in the case of ordinary percolation. The state variable representation also describes the universality 
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class of the percolation models appropriately. 
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Multifractals appear in a wide range of situations like energy dissipation in turbu- 
lent flows 1], electronic eigenstates at metal insulator transition [2], diffusion in porous 
ision limited aggregation J], fluctuations in finance |5(], dynamics of liu- 
6| and many others. The multifractal properties associated with the infinite 
percolation clusters at the percolation threshold Pc is considered in this paper. In ordinary 
percolation (OP) [71], a cluster is generated by occupying a lattice site randomly with a prob- 
ability p or remains empty with a probability (1 — p). Each site of a percolation cluster then 
has two states, occupied or empty. The average q moments of the cluster size distribution 
of percolation clusters are linearly dependent on the moment q and described by a single 
fractal dimension (i/|8|. In order to study the multifractal aspects of percolation clusters, 
usually, a current distribution J9|, [lO| or a random walker ll|] is introduced. However, in the 
presence of rotational constraint on the percolation model, an occupied site has always a 
direction associated with it and a site can be re-occupied from different directions. A site 
can be occupied at most z times from z different directions on a given lattice with coordina- 
tion number z. There are two such well studied rotational percolation models exists, spiral 



percolation SP 12| and directed spiral percolation [13|. In these models, a state variable Sj 
then can be assigned to each lattice site and whose value will correspond to the number of 
times a site is visited during the growth of the cluster. The value of Si then can change 
from to 4 on a square lattice and to 6 on a triangular lattice for SP and DSP models 
whereas it has only two states 1 and in case of OP. In this paper, a new methodology is 
proposed to study the percolation transition in terms of the state variable Sj. Studying the 
physical properties associated with the state variable Sj, the percolation transition is pos- 
sible to establish at the same percolation threshold pc defined geometrically. The spanning 
clusters at p = pc are fractal. Distribution of Si on the fractal substrate is very similar to 
mass distribution on a geometrical support usually taken for multifractal study IJ]. In order 
to explore the multifractal aspects of the spanning percolation clusters in these rotationally 
constrained percolation models aX p = Pc, it is now possible to define a suitable multifractal 
measure in terms of Sj. In this way, one does not need to introduce any other external 
agency like electric current or random walker in the model as it is usually done in the case 
of OP clusters. The variable Sj is inbuilt in the rotational models and represents an inherent 



property of SP and DSP. It is found that the exponents associated with the q moments of 
the measure defined in terms of Si are not hmited by any hnear dependence on the moment 
q for the rotational percolation models. It then indicates that the measure has multifractal 
character. 

Below, description of the rotational percolation models will be given and the multifractal 
aspects of the spanning clusters will be investigated. 

II. ROTATIONAL PERCOLATION MODELS 



There are two well studied rotationally constrained percolation models, spiral percolation 



(SP)[l2l and directed spiral percolation (DSP)[13|]. In these models, clusters are grown 
following single cluster growth Monte Carlo (MC) algorithms. In these algorithms, the 
central site of the lattice is occupied with unit probability. The nearest neighbors of the 
central site is occupied with equal probability p in the first time step. As soon as a site is 
occupied, the direction from which it was occupied is assigned to it. Lists of eligible sites for 
occupation in the next MC time steps are identified and they are occupied with probability 
p. In these algorithms, an occupied site can be reoccupied from a different direction but it is 
forbidden for occupation from the same direction. Once a site is rejected for occupation it is 
forbidden for occupation throughout the simulation from any possible direction. Each MC 
time step can be considered as parallel update of nearest neighbours of already occupied 
sites. The growth of a cluster stops if there is no eligible site available for occupation. 
Identification of eligible sites for occupation in a MC step for SP and DSP are given below. 
In SP, only a rotational field B is applied perpendicular to the plane {xy) of the lattice and 
directed along the negative z-axis (of a right handed coordinate system) . As an effect of the 
B field, empty sites in the forward direction and in the clockwisely rotational direction are 
eligible for occupation. The forward direction is the direction from which the present site is 
occupied and the sense of rotational direction is defined with respect to the forward direction. 
The selection of eligible sites for occupation in SP model is demonstrated in FigHJ^a) for 
both the square and triangular lattices. There are two eligible sites on the square lattice 
and three eligible sites on the triangular lattice due to higher coordination number. The 
eligible sites are then occupied with probability p and the clusters grow isotropically on the 



lattice 



12|. 



In case of DSP, a crossed directional field E is also applied in addition to the rotational 
field B. The E field is applied from left to right in the plane of the lattice and B is, as usual, 
directed along the negative z axis. Due to E field, empty site on the right of an occupied 
site is eligible for occupation along with the eligible sites for occupation due to B field. The 
selection of eligible sites for occupation in DSP model is demonstrated in Fig{T](6) for both 
the square and triangular lattices. The components of the directional constraint remain the 
same on both the lattices whereas there is an extra component of rotational constraint on 
the triangular lattice as in SP. The eligible sites are then occupied with probability p. Once a 
site is rejected for occupation it is forbidden for occupation throughout the simulation from 
any possible direction due to both E and B fields. Because of the simultaneous presence 
of both the E and B fields crossed to each other, a Hall field Eh appears in the system 
perpendicular to both E and 5. As a result, an effective directional constraint E'eff acts on 
the system along the diagonal from left upper to right lower corner of the lattice. Here, the 



clusters grow anisotropically along the effective field -Eeff 13|. 

The coordinate of an occupied site in a cluster is denoted by (x,?/). Periodic boundary 
conditions are applied in both directions and the coordinates of the occupied sites are ad- 
justed accordingly whenever the boundary is crossed. At each time step the span of the 
cluster in the x and y directions L^ = Xmax — Xmin and Ly = Umax — Umin are determined. If 
Lx or Ly > L, the system size, then the cluster is considered to be a spanning cluster. 

Since an occupied site can be reoccupied from different directions due to the presence of 
rotational constraint in both the models, a site then can be occupied at most z times from 
z possible directions on a given lattice where z is the coordination number of the lattice. 
The value of 2; is 4 on the square lattice and it is 6 on the triangular lattice. This is unlike 
in the case of the ordinary percolation where a site is occupied only once. It is an essential 
and a special feature of SP and DSP models. It is thus possible to assign a state variable 
Si to each site and the value of s, will represent the number of times a site is occupied or 
visited. Initially the values of Sj are all set to zero. As soon as a site is occupied from any 
direction, the value of Sj is increased by 1. The values of the state variable is then given by 
Si = 0, ■■■ ,z on a. given lattice. Sj = corresponds to unoccupied site. The state variable 
Si here is similar to the Ising spins with {z + 1) states. In case of OP, Si could have values 
only and 1 corresponding to unoccupied and occupied sites. 

Percolation clusters are generated here in the presence of rotational constraint on the 



square and triangular lattices of size L = 1024 for both SP and DSP models. The per- 
colation transition will be established first at the original pc by calculating "spontaneous 
magnetization" in terms of Sj. The multifractal aspects of the measure distribution on the 
spanning clusters, defined in terms of Si, are then investigated for both isotropic SP and 
anisotropic DSP clusters at p = pc- The cluster properties are averaged over 5 x 10"^ spanning 
clusters. 

III. PERCOLATION THRESHOLD AND SPANNING CLUSTERS 

Geometrically, the critical percolation probability pc is the maximum probability below 
which no spanning cluster appears. At p = Pc, a spanning cluster appears for the first time 
in the system. In single cluster growth approach, the threshold Pc is generally identified by 
measuring the probability to have a spanning cluster (Poo) as a function of p, the occupation 
probability. Pqo goes to zero continuously as p tends to Pc from above. In the state variable 
formalism, the value of pc can be identified in terms of "spontaneous magnetization" M{p) 
defined in terms of the state variable Sj. M{p) is defined as 

. Ntot . L^ 

*°* j=l i=l 

where L is the lattice size and Ntot is the total number of spanning clusters generated. 
mj{p) represents the magnetization per site for the j'th spanning cluster generated using 
single cluster growth method. At p = 1, all the sites of an infinite cluster are expected 
to be occupied z times where z is the coordination number of the lattice and the size of 
the infinite cluster will be of the order of L^, square of the system size. Thus, M(l) is 
equal to z. As p tends to the percolation threshold Pc from above, M{p) is expected to go 
to zero continuously from its maximum value z at p = 1 not only because the sites will 
be occupied less number of times but also the spanning cluster will disappear. M{p) is 
measured on the square lattice {z = 4) for DSP model and it is plotted against p in FiglS It 
can be seen that it is going to zero at p = Pc as expected and the value of Pc is the same as 
that determined by geometrical approach p^ = 0.6550 ± 0.0005 13|. It is also expected that 
M{p) becomes singular aX p = Pc with an exponent (3 as M{p) ^ {p — PcY . In the inset of 
FiglJl the power law has been verified and the exponent (3 is determined approximately as 
(3 ~ 0.32, close to the already obtained value (0.31 ±0.01)|13i]. The state variable formalism 



of the rotationally constrained percolation models is then consistent with that of the usual 
geometrical approach. The value of pc has also been recovered within error bar in the case 
of SP. Other critical properties of the models can also be identified in terms of the state 
variable Si and a scaling theory is possible to develop. 

Typical spanning clusters at p = pc generated on the square lattice of size L = 2^ are 
shown in Figl3I^a) for SP and in Figl3I^6) for DSP. Different values of the state variable Sj 
is represented by different symbols as Sj = (white space), Sj = 1 (filled circle), Sj = 2 
(plus). Si = 3 (filled triangle) and Si = 4 (filled square). It can be seen that not only 
the state variables Sj are randomly distribution of over the spanning cluster but also the 
higher states form small islands allover the spanning cluster. This is similar to the situation 
of mass distribution on a geometrical support generally taken for multifractal study IJ] . 
However, the state distribution over the spanning cluster is not a simple iterative process 
of mass distribution over a geometrical support. The spanning clusters consist of subsets of 
sites occupied once, twice upto a subset of sites occupied z times where z is the coordination 
number of the lattice considered. A particular subset may appear several times in a spanning 
cluster at different stages of growth of the cluster during a large number of MC steps. A 
multiplicative cascade of these subsets is then formed in a complicated manner during the 
growth of the cluster. The Si distribution on the spanning cluster is then expected to have 
many folds. It is then interesting to investigate the moments of the Si distribution over fractal 
objects, the spanning clusters here. It could also be noted that the SP cluster is compact 
and isotropic but the DSP cluster is highly rarefied and anisotropic. The elongation of the 
DSP cluster is along the effective field -Egff appeared in the system. However, the clusters 
are not merely DP clusters along Epff. It has already been found that both SP and DSP 
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belong to new universality classes than that of DP or OP 1121. I13I I. The fractal dimension 

n n 

df of the spanning clusters were found as df ^ 1.733pJ| for DSP and ^ 1.957p^ for SP. 
Geometrical properties of the percolation clusters are governed by this single exponent df. 
However, in the following it will be demonstrated that a measure defined in terms of the 
state variable Si is not restricted by a single exponent rather needs a sequence of fractal 
dimensions to characterize the measure. 



IV. MULIFRACTALITY 



In order to study multifractality a suitable measure has to be defined. In general, a 
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measure is related to the distribution of a physical quantity on a geometrical 



The geometrical support here is the spanning percolation cluster at p = pc for 
the rotationally constrained percolation models. The distribution of the relative probability 
of a state over fractal substrates is a possible multifractal measure here. It is similar to the 
mass distribution on a geometrical support. The multifractal measure /ij is then defined as 

fJ'i = Si/^Si (2) 

where Sj is the state variable associated with each lattice site, /ij can be called as relative 
state variable. Note that the measure /ij is normalized to unity when summed over the 
whole lattice. The maximum value of the measure is fimax = z/ ^Si and the minimum 
non-zero measure is ^mm = VS-^i where z is the coordination number of the lattice. To 
obtain the multifractal nature of the distribution /ij, it is necessary to study the scaling of 
the g- moments of the measure over different length scales on the geometrical support. If 
the measure /ij is multifractal and the support is divided into n^ boxes of size e, then the 
weighted number of boxes N{q,e) is given by 

iV(g,e) = ^/x^^ ;«£--('') (3) 

j=i 

where fij is the sum of the relative state variable in the jth box. Here r(g) can be called as 

"state exponent". The weighted number of box N{q,e) is determined as a function of the 

box size e using box counting method for a given q. The boxes with at least one occupied site 

are only considered. The weighted number of boxes N{q, e) are plotted against the box size e 

for q = —5 to g = +5 for SP in Figj4](a) and for DSP in Fig|4](6) generating spanning clusters 

on the square lattice of size L = 1024. It can be seen that the slopes of the plots change 

continuously for positive g up to g = 0. For g < 0, it seems that the usual box counting 

method adopted here is not working. The values of A^(g, e) remain unchanged over several 

box sizes e starting from the smallest box size for a given g in both the models. It is expected 

that the plot should follow a straight line passing through the points at e = 1 and e = 2^°, 

the system size, in log-log scale since these two extreme points are not effected by the box 

size. It is shown by dashed lines for g = —5 in both the plots. It is observed that the value 



of N{q, e) jumps suddenly when the box size is reduced less than one quarter of the system 
size. This is due to the fact that at this box size there is at least one box appearing with 
a small measure and consequently the sum in Eql3] diverges due to negative moment. The 
appearance of large box sizes with small measures is because of the fact that the spanning 
percolation clusters contain holes of all possible sizes. Difficulties in determining weighted 
number of boxes for g < for the measure distribution on random structures are already 
reported in the literature [l5|. The slopes of the plots in Fig H] also remain almost unchanged 
with the moment q for g < 0. The weighted number of boxes has increased proportionally 
with higher negative moments. It has been verified that the estimation of r(g) by fitting 
only through the smaller box sizes leads to a discontinuity in the plot of r(g) versus q which 
is expected to be a smooth function of q. Discontinuity in the plot of r(g) versus g was 
also observed in the cases of resistance fiuctuations in randomly diluted networks 9| and 



in diffusion limited aggregation 16|. In these cases, there are breakdown of multifractal 
characters for negative moments due to exponential decay of the smallest measures. 

Multifractal characteristics of the spanning clusters of rotationally constrained percola- 
tion models are then analyzed here taking large positive moments, changing q from to 
32. Analysis has been made on the square and triangular lattices of size L = 1024 for both 
SP and DSP models and the results are compared with that of the OP model. In FigJHl 
r(g) is plotted against q, (a) for SP model and (fe)for DSP model. In both the plots the 
squares represent the square lattice data and the triangles represent the triangular lattice 
data. Circles represent the data obtained for OP model. It is found that r(0) is ~ df, the 
fractal dimension of the corresponding spanning clusters and r(l) is ~ here for all three 
models, OP, SP and DSP. r(0) corresponds to the dimension of the support which are the 
spanning percolation clusters of different models considered here and r(l) is zero because 
'^ifj'i = 1. It is interesting to notice that the values of r(g) for DSP and SP model depend 
on the moment g in a nonlinear way for positive moments. If the measure /ij is characterized 
by a single fractal dimension df, the state exponent r(g) should have a constant gap between 
two consecutive exponents|8| and consequently should have a linear dependence on g. In 
that case, a relationship between r(g) and g in terms of df can be obtained as 

riq) = -iq-l)df. (4) 

This relation is shown in Fig|5] by a solid line for OP taking df = 91/48 ~ 1.896. The 



values of r(g) obtained numerically for OP (circles) considering a two state model for which 
Si = and 1 are in good agreement with with EqJH There are few more things to notice. 
First, the state exponents of OP obey the constant gap equation given in EqJUas expected. 
The constant gap scaling was also observed for the mean number of distinct sites visited by 



a random walker on spanning percolation cluster by Murthy et a/[l7j. As a consequence, 
the measure distribution is mono-fractal. Second, the values of r(g) for SP and DSP are 
deviated form straight line behaviour and have non-linear dependence on the moment q. 
Thus, each moment of the measure /ij needs a new exponent to characterize in these models. 
Third, the functional dependence of r(g) on q is found different for all three models. This 
is expected because the models, OP, SP and DSP, belong to different universality classes. 
Fourth, the r(g) values on the square and triangular lattices are almost the same for the SP 
model whereas they are considerably different for the DSP model. This is also in agreement 
with the fact that the critical properties hold universality in the SP model whereas there is 
a breakdown of universality in DSP model between the square and triangular lattices in two 
dimensions 19|]. Finally, the fact that a sequence of exponents is required to characterize the 
moments of the measure confirms the multifractal nature of /ij distribution in SP and DSP 
models. It should be noted here that in SP and DSP the spanning clusters consist of four or 
six subsets depending on the number of nearest neighbours on a given lattice. The values 
of r(g) then may be possible to obtain in terms of the fractal dimensions of the subsets 
consisting the spanning cluster coupled with a nonlinear dependence on q. However, it is 
difficult to determine the fractal dimensions of the individual subsets as they are generally 
small isolated islands in a spanning cluster as well as the nonlinear nature of T{q). 

The associated fractal dimensions /(a) with the measure and the corresponding Lipschitz- 



Holder exponent a can be obtained through a Legendre transformation 18| of the sequence 
r(g). The Legendre transformation is given below 

«(g) = -^, /(«) = g«(g) + r(g). (5) 

The fractal dimensions /(a) are plotted against a in Fig|6l The values of /(a) obtained 
for SP clusters are shown in Figl6] (a) and that of the DSP clusters are shown in Figl6] (b). 
Since in the case of OP, the state exponent follows a constant gap equation (Eq|4]) it is 
expected that /(«) versus a will be represented by a point f{a) = a = dj. It is shown by 
an open circle in FiglHl It has been verified measuring slopes at different regions of r(g) 

9 



versus q for OP that the slopes remain within the error bar of the point shown. In the 
cases of SP and DSP, spectra of /(a) against a are obtained because r(g) has a non hnear 
dependence on q. The symbols square and triangle correspond to the square lattice and 
the triangular lattice data respectively. There are few things to observe. First, fmax{ct) 
corresponds to df of the respective models. Second, the /(«) curves are always < df since 
the supports are spanning percolation clusters of fractal dimension df. Third, the spectrum 
of fractal dimensions /(a) for SP and DSP are found different. It means that not only 
the mass fractal dimension df is different but also the whole set fractal dimensions /(«)s 
are different. It is expected because SP and DSP belong to different universality classes. 
Fourth, in case of SP, the spectrum of /(a) for the square lattice is identical with that of the 
triangular lattice for lower moments (starting from the same mass fractal dimension) and 
slightly different at large positive moments. In case of DSP, the spectra of /(a) on the two 
lattices are considerably different over the full range of positive moments considered here, 
starting from two different mass fractal dimensions. This again confirms the universality of 
critical exponents in SP and breakdown of universality of critical exponents in DSP between 
square and triangular lattices in two dimensions [1 91]. Fifth, the values of /(amj„) for both SP 
and DSP clusters are not equal to zero. This means that in these cases, the rarest regions 
of measure fimax distribution are still fractal in the limit q -^ oo. It is evident from the 
spanning clusters configurations given in FigOJthat the fimax distribution appears in small 
islands allover the spanning clusters and not as a point distribution. The fractal dimension 
of fimax distribution of SP clusters is found little higher than that of DSP clusters. This 
is due to the presence of an extra directional constraint in the DSP model which takes the 
growth of the cluster away from a fimax point whereas due to pure rotational constraint in 
the SP model the probability of growth around a fimax point is little higher in comparison 
to DSP. It can also be noticed that, in case of DSP, the whole /(«) spectrum is shifted 
upward in order to match with the mass fractal dimension of the spanning DSP clusters on 
the triangular lattice. However this is not the case in SP model. Finally, the values of /(a) 
converges at a minimum of Lipschitz- Holder exponent amm- The Lipschitz-Holder exponent 
a(^) is defined as fx^ = S"^^^ where /i^ = fi{^ + 5) — ii{^) is the increment in the measure 
over a length ^ to ^ + 5 [iJ]. The amin value corresponds to the minimum of ^, the length 
scale associated with fimax clusters in this case. It could also be noticed that in both SP 
and DSP, the values of amm on the triangular lattice is found little higher than that of the 
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square lattice. This is due to the fact that the number of fimax points are generally higher 
on the triangular lattice than that on the square lattice. 

V. CONCLUSION 

Using the concept of state variable, the percolation transition in rotationally constrained 
models is established at the same percolation threshold pc defined geometrically. A relative 
state variable is defined to study the multifractal aspects of the spanning clusters at p = Pc- 
It is found that the g-moments of the measure is characterized by a sequence of "state 
exponents" r(g) for both SP and DSP. The existence of a sequence of state exponents 
confirms the multifractal character of the distribution of relative state variable on the infinite 
clusters of SP and DSP. The OP spanning clusters are not found multifractal in this measure. 
Taking Legendre transformation of r(g), different spectra of associated fractal dimensions 
f{a) as a function of Lipschitz-Holder exponents a are obtained. The universality of critical 
exponents in SP and breakdown of universality in DSP are also confirmed by the multifractal 
spectrum of fractal dimensions. The formalism of state variable is thus found suitable for 
studying percolation transition and multifractal aspects of certain percolation models. 
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FIG. 1: Selection of empty nearest neighbours in a MC step in (a) SP model and (6) DSP model 
on the square and triangular lattices. Black circles are the occupied sites and the open circles are 
the empty sites. Thick long arrows from left to right represent directional constraint (E). The 
clockwise rotational constraint (B) is shown by encircled dots. The central site here is occupied 
from site 2 and shown by short thick arrows. The eligible sites for occupation due to E field are 
shown by dotted arrows and thin solid arrows indicate the same due to B field on both the lattices. 



13 



2.5 
2.0 
S 1-5 
1.0 
0.5 



0.0 




>t 



0.65 



6 , -5, ,-4 



0.70 



0.75 
P 



0.80 



FIG. 2: Plot of spontaneous magnetization M{p) against p for DSP model defined on a square 
lattice of size L = 1024. Percolation threshold pc is marked by a cross on the p axis. In the inset, 
M is plotted with {p — Pc)- From the slope, the exponent /3 is obtained as w 0.32. 
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FIG. 3: Typical spanning clusters at p = pc on the square lattice of size L = 2^ (a) for SP and 
(b) for DSP models. The encircled dots represent the rotational field B and the arrows represent 
the directional field E. Different symbols in the clusters represent different values of Si as (a) for 
Si = 1, (+) for Si = 2, (a) for Sj = 3 and (■) for Si = 4. The empty white space represents Sj = 0. 
It can be seen that the state variable is randomly distributed over the fractal spanning clusters. 
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FIG. 4: Plot of weighted number of box N{q,e) versus the box size e for q = —5 to g = 5 in step 
of 1 for SP in (a) and for DSP in (b) for the spanning clusters generated on the square lattice of 
size L = 1024. The symbols are: (•) for q = —5, (n) for q = —4, (o) for q = —3, (A) for q = —2, 
(<l) for q = -1, (V) for q = 0, {>) for q = 1, (+) for q = 2, (x) for q = 3, (*) for g = 4, and (o) for 
q = 5. It can be seen that box counting method is not working for g < 0. The expected behaviour 
is shown by dashed lines for g = — 5 in both (a) and (6). 
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FIG. 5: Plot of the state exponent r(g) versus the moment q for SP in (a) and for DSP in (b) 
for q > 0. Squares represent the square lattice data and triangles represent triangular lattice data 
respectively. Circles represent data of OP spanning clusters on a square lattice. The solid straight 
line represents the linear dependence of r((/) on q expressed in Eq|3]for OP. The measured values of 
T(g) for OP follows the straight line behaviour. For DSP and SP, T{q) has a non-linear dependence 
on the moment q. 



17 







FIG. 6: Plot of fractal dimension /(a) against the Lipschitz-Holder exponent a for SP (a) and 
DSP (b). Squares represent the square lattice data and triangle represent triangular lattice data 
respectively. For OP, it is a single point at f{a) = a = dj and represented by a circle. For DSP 
and SP, different spectra of f{a)s are obtained. The spectra on the square and triangle lattice 
differ considerably for DSP whereas for SP they are almost identical. 
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